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ABSTRACT 

We use the Riess et al. (2004) supernova la apparent magnitude versus red- 
shift data and the Allen et al. (2004) galaxy cluster gas mass fraction versus 
redshift data to constrain dark energy models. These data provide complemen- 
tary constraints that when combined together significantly restrict model param- 
eters and favor slowly-evolving dark energy density models, close to the Einstein 
cosmological constant limit of dark energy. 

Subject headings: cosmology: cosmological parameters — cosmology: observations — 
X-rays: galaxies: clusters — supernovae: general 



1. Introduction 

Less than a decade ago, astronomers discovered that the expansion of the Universe 
was speeding up (Riess et al. 1998; Perlmutter et al. 1999). In the context of Einstein's 
general relativity, this acceleration is attributed to dark energy that varies slowly with time 
and space, if at all. More recent supernova la (SNIa) redshift-magnitude data,^ and cosmic 
microwave background (CMB) anisotropy data — that indicates the Universe is spatially 



^Present address: Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822; gchen@ifa.hawa- 
ii.edu 

^For the SNIa data see, e.g., Riess et al. (2004, hereafter R04), Nobih et al. (2005), Clocchiatti et al. 
(2006), and Astier et al. (2006). SNIa data constraints on cosmological parameters are also discussed by 
Wang & Tegmark (2005), Lazkoz et al. (2005), Hannestad (2005), Shafieloo et al. (2006), Zhang & Wu 
(2005), Godlowski & Szydlowski (2005), Sethi et al. (2005), Bento et al. (2005), and Jassal et al. (2006). 



- 2 - 



flat (see, e.g., Podariu et al. 2001b; Durrer et al. 2003; Melchiorri & Odman 2003; Bennett et 
al. 2003; Page et al. 2003) — combined with a low mean value for the nonrelativistic matter 
density parameter Qm (Chen & Ratra 2003b, and references therein), strengthens the claim 
that a previously undetected dark energy accounts for a major part of the cosmic energy 
budget. For reviews of the current situation see Peebles & Ratra (2003), Steinhardt (2003), 
Carroll (2004), Padmanabhan (2005), and Perivolaropoulos (2006). 

While CMB anisotropy and SNIa constraints on dark energy will continue to improve,^ 
it is important to use many different tests to constrain dark energy models. This will allow 
for consistency checks, as well as provide complementary constraints on dark energy. A 
number of other tests have been considered. These include the redshift-angular size test 
(see, e.g., Chen & Ratra 2003a; Podariu et al. 2003; Jackson 2004; Puetzfeld et al. 2005; 
Daly & Djorgovski 2005), the strong gravitational lensing test (see, e.g., Fukugita et al. 1990; 
Turner 1990; Ratra & Quillen 1992; Chae et al. 2004; Sereno 2005; Alcaniz et al. 2005), the 
galaxy cluster gas mass fraction versus redshift test (see, e.g., Sasaki 1996; Pen 1997; Allen 
et al. 2004, hereafter A04; Chen & Ratra 2004, hereafter CR; Rapetti et al. 2005; Alcaniz & 
Zhu 2005; Kravtsov et al. 2005; Chang et al. 2006), and the baryon acoustic oscillation test 
(see, e.g., Glazebrook & Blake 2005; Angulo et al. 2005; Wang 2006). In addition, structure 
formation in dark energy models will soon provide a useful diagnostic of model parameters 
(see, e.g., Horellou & Berge 2005; McDonald et al. 2006; Le Delliou 2005; Bartelmann et al. 
2005; Percival 2005; Mainini 2005). In this paper we constrain dark energy models by doing 
a joint analysis of the SNIa apparent magnitude versus redshift data of R04 and the galaxy 
cluster gas mass fraction versus redshift data of A04. 

There are a number of different dark energy models under current discussion.^ The cur- 
rent "standard" cosmological model is ACDM (Peebles 1984), where the energy density of 
the low-redshift Universe is dominated by a cosmological constant A (with a constant energy 
density p\), with nonrelativistic matter — mostly cold dark matter (CDM) — playing a 
subdominant role. Models in which the dark energy is due to a slowly varying, in time and 
space, scalar field (0), have also attracted recent attention. A prototypical 0CDM example 



^For near future SNIa constraints sec SoUerman et al. (2005) and Hamuy et al. (2006). Within 
a decade or so, the proposed JDEM/SNAP space mission should provide even tighter constraints (see 
http://snap.lbl.gov/ and Podariu et al. 2001a; Crotts et al. 2005; Albert et al. 2005, and references therein). 
Mukherjee et al. (2003), Caldwell & Doran (2004), and Lee et al. (2006) discuss CMB constraints on dark 
energy. 

^Recent papers include Ncupanc & Wiltshire (2005), Pietroni (2005), Zhao et al. (2005), CapozzicUo et 
al. (2005), Franga (2005), Barenboim et al. (2005), Barnes et al. (2005), Abdel-Rahman & Riad (2005), 
Barger et al. (2005), Sola & Stefancic (2006), and Arbey (2006). 
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has scalar field potential energy density V{(f)) oc a > 0, at low redshift (Peebles & 
Ratra 1988; Ratra & Peebles 1988); again CDM plays a subdominant role. Also discussed 
is the XCDM parametrization for time-varying dark energy. This parametrization approxi- 
mates dark energy by a fluid with a negative time-independent equation of state parameter 
w = P/p, where P is the fluid pressure and p the energy density. This is an inaccurate 
approximation during the scalar field dominated epoch when w is time dependent (see, e.g., 
Ratra 1991). We consider spatially flat spacetimes for the 0CDM model and the XCDM 
parametrization, but allow spatial curvature to be a free parameter in the ACDM case. In 
this paper we constrain parameters of these three models by using SNIa and galaxy cluster 
data. 

In § 2 we show how we use the supernova and galaxy cluster data to constrain cosmo- 
logical parameters. Results are presented and discussed in § 3. We conclude in § 4. 

2. Computation 

A04 list x-ray gas mass fractions for 26 rich clusters, determined from Chandra obser- 
vations. These clusters lie at redshifts between 0.08 and 0.89. Following A04 and CR, we use 
this data to determine the probabihty distribution function (likelihood) L'-'{flM,P, h). Here 
h is the Hubble constant in units of 100 kms~^Mpc~^, and p is the cosmological constant 
density parameter Qa for the ACDM model, w for the XCDM parametrization, and a for 
the 0CDM model. To derive L^{Qm,P, h) we marginalize over the bias factor b, as well as 
over fljjh^, where Qi, is the baryonic mass density parameter. Later we will need to marginal- 
ize over h. We account for uncertainties in b, Qi,h'^, and h by using Gaussian priors with 
b = 0.824 ± 0.089, n^h^ = 0.0214 ± 0.002, and h = 0.72 ± 0.08, all one standard deviation 
errors.^ See A04 and CR for more detailed discussions of the procedure we use. 

We also use the SNIa apparent magnitude versus redshift data from R04, in particular, 
the gold data set of 156 SNIa with redshifts up to almost 1.8. From this data we determine 



'''See Ganga et al. (1997) for a discussion of this method. As noted in CR, there is a more precise estimate 
of h than the value quoted above, h = 0.68 ± 0.04 (Gott et al. 2001; Chen et al. 2003, and we have halved 
the quoted two standard deviation errors). Also, the Ofo/i^ value used above is consistent with the estimate 
from the WMAP CMB anisotropy measurement and the primordial deuterium abundance measurement, but 
significantly higher than an estimate based on the primordial helium and lithium abundance measurements 
(see, e.g., Peebles & Ratra 2003; Steigman 2005; Fields & Sarkar 2006), so we also consider the Peebles & 
Ratra (2003) summary estimate fife/i^ = 0.014 ± 0.004. Constraints derived using the values of h and Qfe/i^ 
quoted in this footnote are shown as darker dashed lines in Figs. 2 — 4, while constraints from the priors 
quoted in the main text are shown using lighter continuous lines. 
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the probability distribution function L^{Q,M,P,h), following, e.g., R04 or Podariu & Ratra 
(2000). 

The joint likelihood for the SNIa and galaxy cluster data is the product of the two 
individual likelihoods. The two-dimensional probability distribution function for Q^m and p, 
L(fiM,p), is determined by marginalizing this product over h. For each of the three models 
mentioned above, we compute L(Qm,p) on the two-dimensional {VLm^p) grid. The 1, 2, and 
3 a confidence contours are the set of points where the likelihood is e~^'^°/^, e~^'^^/^, and 
g-ii.8/2 q£ ^YiQ maximum likelihood value. 



3. Results and Discussion 

Figure 1 shows the R04 gold SNIa constraints on the 0CDM model with V{(t)) oc 
determined by marginalizing L^{QM,Ci,h) over h. These contours are significantly more 
constraining than those derived by Podariu & Ratra (2000, Fig. 1) from earlier SNIa data. 

Figure 2 shows the combined SNIa and galaxy cluster constraints on the ACDM model. 
These joint constraints are significantly tighter than those derived using either the SNIa 
data (R04, Fig. 8) or the galaxy cluster data (A04, Figs. 4 and 8; CR, Fig. 1) alone. This 
is because the galaxy cluster data tend to tightly constrain Qm (A04, CR), while the SNIa 
data tend to tightly constrain a linear combination of Qa and Qm (R04). Together, the data 
focus attention on a small part of parameter space near Qm ~ 0.3 and fl\ ^ 0.7 where the 
Universe is spatially flat. 

Figure 3 shows the joint data constraints on the XCDM parametrization. Models close 
to Qm ~ 0.3 and w ^ —1 are favored. We emphasize that w = —1 corresponds to the 
spatially-flat ACDM model. 

Figure 4 shows the joint data constraints on the 0CDM model with V{(f)) oc Models 
close to Qjvf ~ 0.25 and a ~ are favored. We note that a = corresponds to the spatially- 
flat ACDM model. 

The SNIa and galaxy cluster data together favor a spatially-flat ACDM model with 
Qm ~ 0.25 — 0.3. It might be signiflcant that other data also favor this value of Qm (see, 
e.g., Chen & Ratra 2003b; Spergel et al. 2003). We emphasize, however, that slowly-evolving 
dark energy models are also consistent with this data. 
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4. Conclusion 

We use supernova la apparent magnitude and x-ray cluster gas mass fraction data to 
jointly constrain cosmological parameters. These data sets provide complementary con- 
straints and together they tightly constrain cosmological parameters. The spatially-fiat 
ACDM model is favored, although a 0CDM model with slowly-decreasing dark energy can- 
not yet be ruled out. 

More and better SNIa and galaxy cluster data will allow for tighter constraints on model 
parameters. Equally important would be a resolution of the conflict between the two different 
sets of estimates of Qb, as well as tighter hmits on h. All these are hkely to happen in the 
next few years. 

We acknowledge useful discussion with M. Sayler, K. Teramura, and J. Whitmer, and 
support from DOE grants DE-FG03-99ER41093 and DE-FG02-00ER45824, NASA grants 
AISR NAG5-11996 and ATP NAG5-12101, and NSF grants AST-0206243, AST-0434413, 
and ITR 1120201-128440. 
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Fig. 1. — Probability distribution function confidence contours (1, 2, and 3 a from inside 
to outside) for the 0CDM model using the R04 gold SNIa sample. The dot denotes the 
maximum likelihood and is at = 0.30 and a = 0.0. 
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Fig. 2. — Probability distribution function confidence contours (1, 2, and 3 a from inside 
to outside) for the ACDM model. The solid gray lines are computed using h = 0.72 ± 0.08 
and Vlbh^ = 0.0214 ± 0.002, with maximum likehhood at = 0.29 and = 0.75. The 
black dotted lines use h = 0.68 ± 0.04 and Qbh"^ = 0.014 ± 0.004, with maximum likelihood 
at ^Im = 0.23 and = 0.65. The slanting dash-dot line indicates spatially-flat models and 
models to the left and above the sohd black line in the left top corner do not have a big 
bang. 
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Fig. 3. — Probability distribution function confidence contours (1, 2, and 3 a from inside to 
outside) for the XCDM parametrization. The solid gray lines are computed using h = 0.72 ± 
0.08 and Qbh^ = 0.0214±0.002, with maximum likelihood at = 0.29 and w = -1.04. The 
black dotted lines use h = 0.68 ±0.04 and Qbh"^ = 0.014 ±0.004, with maximum likelihood at 
flM = 0.23 and w = —0.88. The horizontal dash-dot line at w = —1 indicates spatially-fiat 
ACDM models. 
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Fig. 4. — Probability distribution function confidence contours (1, 2, and 3 a from inside 
to outside) for the 0CDM model with V{(f)) oc 0^". The solid gray lines are computed 
using h = 0.72 ± 0.08 and n^h^ = 0.0214 ± 0.002, with maximum likelihood at VIm = 0.28 
and a = 0.0. The black dotted lines use h = 0.68 ± 0.04 and ^bh^ = 0.014 ± 0.004, with 
maximum likelihood at Qm = 0.22 and a = 0.45. The horizontal a = axis corresponds to 
spatially-flat ACDM models. 



